Radio Quiet AGNs as Possible Sources of Ultra-high Energy Cosmic Rays 
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Active galactic nuclei (AGNs) have been one of the most widely discussed sources of ultrahigh- 
energy cosmic rays (UHECRs). The recent results of Pierre Auger observatory (PAO) have indicated 
a possible composition change of UHECRs above ~ 10 18 ' 5 eV towards heavy nuclei. We show here 
that if indeed UHECRs are largely heavy nuclei, then nearby radio quiet AGNs can also be viable 
sources of UHECRs. We derive constraints on the acceleration sites which enable acceleration of 
UHECRs to 10 20 eV without suffering losses. We show that the acceleration of UHECRs and the 
survival of energetic heavy nuclei are possible in the parsec scale weak jets that are typically observed 
in these objects, the main energy loss channel being photodisintegration. On this scale, energy 
dissipation by shock waves resulting from interactions inside a jet or of the jet with surrounding 
material are expected, which may accelerate the particles up to very high energies. We discuss the 
possible contribution of radio-quiet AGNs to the observed UHECR flux, and show that the required 
energy production rate in UHECRs by a single object could be as low as « 3 x 10 39 ergs -1 , which 
is less than a percent of the bolometric luminosity, and thus energetically consistent. We discuss 
consequences of this model, the main one being the difficulty in detecting energetic secondaries 
(7-rays and neutrinos) from the same sources. 

PACS numbers: 98.54.Cm, 95.85.Ry, 98.70.Sa 



I. INTRODUCTION 

The origin of ultrahigh-energy cosmic rays (UHECRs) , 
cosmic rays with energies above > 10 18 5 eV, is still a 
mystery. It is commonly believed that the sources of 
UHECR's are extragalactic, since at these energies the 
cosmic rays cannot be confined by the magnetic field of 
our galaxy (for recent reviews, see Refs. 0, Q). Due to 
energy losses by photo-meson production with the cosmic 
microwave background (CMB), the sources of UHECR's 
are further limited to a distance < 100 Mpc from the 
observer (this limit is known as the 'Greisen-Zatsepin- 
Kuzmin (GZK) cutoff; Refs. 0, Q). While early ob- 
servations [5( suggested a violation of the GZK cutoff at 
high energies, recent, high statistics observations by the 
High-Resolution Fly's eye (HiRes) experiment and the 
Pierre Auger observatory (PAO) finds a downturn in the 
spectrum consistent with the GZK predictions [1, 0] • 

Extracting information from observations is, however, 
not easy. For example, observations of the spectrum it- 
self did not lead, so far, to a firm conclusion about the 
origin of the steepening at high energies (i.e., whether 
it originates from the GZK cutoff or not), and therefore 
other possibilities are still allowed. 

One of the consequences of this difficulty is the lack of 
a clear theoretical picture of the extragalactic UHECR 
sources. Possible candidates of UHECR sources are lim- 
ited among the known astrophysical objects, because 
of the difficulty in satisfying the two nearly contra- 
dictory requirements: A strong magnetic field in the 



source is needed to confine the accelerated cosmic rays, 
while the magnetic and photon fields cannot be too 
strong in order to avoid too much synchrotron radia- 
tion and photo- meson energy losses [||. Several possi- 
ble sources of UHECRs that fulfill the constraints men- 
tioned above are often discussed in the literature. The 
most widely discussed candidates are gamma ray bursts 
(GRB )[|, Gl, [p, p. 111 03, low luminosity GRBs 
and hypernova e 114. \WL flol . and jets in radio-loud (RL) 
AGNs [13, EE [TaHlallJla HI- UHECR acceleration 
in the vicinity of black holes has also been considered 
for AGNs, including radio-quiet (RQ) AGNs (see, e.g., 
Refs. [241. l25ll ). Additional suggestions include magne- 
tars dalljand clusters of galaxies [HIH- 

An important step towards understanding the origin of 
UHECRs was taken recently with the discovery of a cor- 
relation between the arrival direction of the highest en- 
ergy cosmic rays and nearby galaxies in the [30] (VCV) 
catalog of active galactic nuclei (AGNs) [12th edition] 
[3lL I32I ]. These results, however, should be taken with 
great care. It was noted already by Refs. [3l],[32|] that the 
observed correlation alone is insufficient to conclude that 
the sources of UHECRs are AGNs, because AGNs them- 
selves are clustered within galaxies, and thus may only 
be the tracers of the true sources. Indeed, a correlation 
was found between the arrival directions of UHECRs and 
galaxies [33|,[34],[35|], which suggests a correlation with the 
large scale structure (LSS) in the nearby universe. 

The emerging picture gets further complicated by the 
recent results of HiRes experiment which do not confirm 
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the correlation found by PAO [3(|. Moreover, the latest 
PAO results [37| suggest that the correlation degree is 
reduced compared to that obtained by earlier measure- 
ments. 

An additional clue on the nature of UHECRs sources 
may come from measurements of their chemical compo- 
sition. However, deducing the chemical composition of 
UHECRs from current observations is also difficult due 
to our poor knowledge on hadronic interactions. PAO re- 
sults indicate a heavy composition at the highest energies 
[H, Hs| , while HiRes results suggest a proton composi- 
tion (see Ref. Q, and references therein). The results 
of these two experiments are therefore clearly in contra- 
diction. Moreover, the heavy element composition sug- 
gested by PAO results may be in contradiction with the 
anisotropy result of the same experiment (see below) . 

The chemical composition of UHECRs may be crucial 
in determining the nature of their source. If UHECRs 
are protons, the theoretical restrictions on the physical 
conditions inside potential sources are rather extreme, so 
that if indeed AGNs are the sources of UHECRs, then 
only powerful AGNs such as RL AGNs can be viable 
sources. In addition, the known values of the galactic 
magnetic field, as well as observational upper limits on 
the intergalactic magnetic fields (IGMFs) in voids [i(| 
and theoretical investigations about IGMFs in the struc- 
tured region [4l[ allow one to expect a correlation be- 
tween the arrival directions of UHECRs and AGNs [13] , 
which is tentatively suggested by PAO. On the other 
hand, a drawback of the idea that AGNs are sources of 
UHECRs is the fact that a strong correlation with RL 
AGNs is not found, except for very few cases, such as 
Cen A or Cen B [43|. In addition, a correlation with 
the most powerful AGNs, FR II galaxies, is currently not 
confirmed. Moreover, the number density of FR II galax- 
ies seems too small to explain the observed small scale 
anisotropy Furthermore, the observed bolometric 

luminosity of AGNs which are correlated with the arrival 
direction of UHECRs seems too small to satisfy the min- 
imum conditions for UHECRs acceleration in continuous 
jets [H| . These results indicate that only very few nearby 
AGNs are sufficiently powerful to accelerate particles to 
the observed ultrahigh energies. 

In order to overcome this problem, it was suggested by 
Ref. [46| that UHECR production may occur transiently 
during giant AGN flares caused by tidal disruptions. Al- 
though this may be a valid scenario, several restrictions 
on the duration and rate of such flares can be derived 
[47l |48| . Since the actual activity of AGNs is not known, 
concrete conclusions cannot be drawn at this stage (see 
also Ref. (H). 

An alternative picture emerges if UHECRs are com- 
posed of heavier elements, such as iron nuclei. The am- 
biguity in the interpretation of the data led [1(| HU to 
study a mixed composition scenario above - I0 185 eV. 
A similar model (albeit with somewhat lower maximum 
energy of UHECR) was studied recently by [12] • In these 
models, significant correlation between the arrival direc- 



tions of UHECRs and galaxies is not expected, because 
of the strong deflection of heavy nuclei in galactic mag- 
netic field. In addition, protons dominate the cosmic 
rays spectra only up to energies ~ I0 18 ' 5 eV, which im- 
plies that the requirement on the physical conditions at 
the acceleration site of potential sources is significantly 
loosened. 

As we show here, a heavy element composition of UHE- 
CRs opens the possibility for RQ AGNs, which compose 
the majority of the AGN population within ~ 100 Mpc, 
to be the sources of UHECRs. As opposed to RL AGNs, 
in which evidence exist for strong jets on > kpc scale, 
in RQ AGNs there are no evidence of such strong jets 
[i)2l ]. Nonetheless, in recent years there has been an accu- 
mulation of evidence for weak, parsec-scale jets in these 
objects (Refs. [13, [H, HE HE HI, HI, and references 
therein). Internal interactions of blobs inside a jet, or 
of the jets with their environment, produce shock waves, 
which can accelerate particles. Alternatively, cosmic rays 
may be accelerated by the dissipation of magnetic energy 
inside these jets (see, e.g., Ref. |60J ) . Moreover, as was 
pointed out by Ref. (6fJ , particle acceleration could take 
place in shocks within the corona, which could be caused 
by abortive jets. 

In this work we thus discuss weak jets in RQ AGNs 
as possible sources of UHECRs, under the assumption 
that at very high energies the composition of UHECRs is 
dominated by heavy (possibly, but not necessarily, iron) 
nuclei. As we will show below, RQ AGNs can fulfill the 
requirements needed from sources of UHECRs under this 
assumption. We assume here that particles acceleration 
takes place inside the jets, so that this work is different 
from previous works that considered UHECR production 
in the vicinity of black holes that do not have jets (e.g., 
Ref. [H). 

This paper is organized as follows. In fjll] we present 
the basic constraints on the physical conditions at the 
acceleration sites of heavy nuclei. In particular, in mi Al 
we derive the constraints from photodisintegration and 
photomeson production. We show that these conditions 
can be fulfilled easily. In mil] we discuss the efficiency 
of particle acceleration required to explain the observed 
flux of UHECRs on earth, under the assumption that RQ 
AGNs are the only sources of UHECRs. We further de- 
rive constraints on the value of the magnetic field at the 
acceleration site, and show that a lower limit of few per- 
cents of equipartition value is required. We summarize 
and conclude in $V] 



II. CONSTRAINTS ON THE EMISSION SITE 
OF UHECRS 

The acceleration sites of UHECRs are required to fulfill 
two basic conditions (see, e.g., Refs. d,[l3,l3|)- The first 
condition is that the accelerated particles should be con- 
fined to the acceleration region (this is also known as the 
"Hillas condition", see Ref. [8|). Assuming that the ac- 
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celeration results from electromagnetic processes within 
an expanding plasma, this condition is equivalent to the 
requirement that the acceleration time is shorter than the 
dynamical time scale. Considering jetted plasma, which 
moves relativistically with velocity (3c, or Lorentz factor 
r = (1 — /3 2 ) -1 / 2 , the acceleration time scale in the co- 
moving frame can be written as i acc ~ i]E ob /(TZqBc) 
[93l | . Here, E ob = TE' is the observed energy of the par- 
ticle [94| . E' is its energy in the comoving frame, B is 
the magnetic field, and Zq is the nucleon charge. The 
factor 77 > 1 is a dimensionless factor, whose exact value 
is determined by the (yet uncertain) details of the accel- 
eration mechanism. For example, in the non-relativistic 
diffusive shock acceleration mechanism, this factor corre- 
sponds to ij — (20/3)(3~ 2 in the Bohm limit for parallel 
shocks (e.g., Refs. [H, H3|). For relativistic shocks, r\ ~ 
a few can be expected but larger values are also possible. 
The condition £ acc < id yn = r/T/3c, gives 
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r)E ob (3 _ 3.3 x 10 17 r]P 
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(1) 



where r is the radial distance from the source at which 
particle acceleration takes place, measured in cm. Here 
and below, we use the convention Q = 10 x Qx in cgs 
units. 

Writing the energy density in the magnetic field as a 
fraction eb of the photon energy density at the accelera- 
tion site ub = B 2 /8tt = eg«, where u = L / (Aitr 2 T 2 (3c) is 
the photon energy density (representative of the electron 
energy density) and L is the photon luminosity, equation 
[T]can be written in the form (Ref. [46j]) 



e B L > 1.7 x 10 45 rfZ- 



T 2 (3 3 E°° - ergs 



oh- 



(2) 



Equations [TJ [2] give basic conditions required from a 
source capable of accelerating particles to the maximum 
observed UHECR energy, ~ 10 20 eV. However, as the de- 
tails of the acceleration process are uncertain, the actual 
restriction may depend on additional phenomena, such 
as the details of particle escape from the acceleration re- 
gion. The details of particle escape are unknown, and 
the maximum energy may be escape-limited if the es- 
cape time scale is shorter than the dynamical time scale. 
In this paper, we simply assume that acceleration is effi- 
cient, such that £ C sc ^ idyn, which can be expected when 
the Bohm limit is achieved over the size of r/T. 

For protons (Z=l), equation [2] necessitates a mag- 
netic luminosity which is 1-2 orders of magnitude above 
the bolometric photon luminosity observed in nearby 
AGNs associated with UHECR's in the VCV catalog, 
L boi ~ 10 44 ergs _1 [45]. We therefore expect that for 
proton composition of UHECR, acceleration to the high- 
est energies is possible only in powerful RL AGNs (for 
which L > 10 44 ergs -1 ). However, for heavier nu- 
clei, equation [2] implies a much weaker constraint: e.g., 
for Carbon nuclei (Z = 6) one obtains e^L > 5 x 
10 43 r/ 2 r 2 /3 3 E ob 2 ergs -i_ p or iron nuclei, Z = 26 
and the constraint on the luminosity further eases to 



e B L > 2.5 x 10 4 Vr 2 /3 d E$f ergs" 1 , well within the 
limits of many nearby AGNs' observed luminosity (ex- 
cept for low-luminosity AGNs, with L< 10 42 ergs -1 ) [9a ]. 

The second condition is that the acceleration time is 
shorter than all the relevant energy loss time scales such 
as synchrotron loss time. For the synchrotron loss time, 
we have tcooi^n = (6Trm 4 c 3 rA 4 )/(a T m 2 e B 2 E ob Z 4 ) (see, 
e.g., Ref. [65j|). Here, ot is Thomson's cross section, 
m p and m e are the proton and electron masses and Am p 
is the mass of the nucleon (for iron nuclei, A = 56 and 
Z = 26). The requirement i acc < i C ool,syn results in 

E ob < 2 x 10 20 T rT 1 / 2 £T 1/2 A 2 Z^' 2 eV. (3) 

For iron nuclei, equation [3] gives E ob < 5 x 
10 21 Tr]- 1 / 2 (B/lG)- 1 / 2 eV, which implies that in order 
to enable acceleration of cosmic rays to the highest ob- 
served energies, ~ 10 20 eV, the strength of the magnetic 
field at the acceleration site should not exceed B < few - 
few tens G (as long as T ~ 1, as is expected in RQ AGNs. 
See Ref. H3). 

The constraints of equations [T] and [3] can be combined 
to constrain the radius of the acceleration site, 

r acc Er>9x 10 16 J] 2 T~ 2 (3 Z 2 A- 4 E°f cm. (4) 

For iron nuclei, equation 0] implies r acc > 6 x 
10 12 if T -2 (3E^ cm, which is comparable to the 
Schwarzshield radius of typical black holes in AGNs 
(with characteristic mass Mbh = 10 8 M Q ), rs c h.M s = 
3 x 10 13 cm. 

Comparison of the synchrotron cooling time and the 
dynamical time shows that if the acceleration takes place 
at radii smaller than 
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(5) 



then the synchrotron cooling time for the most en- 
ergetic particles is shorter than the dynamical time, 
icooi,syn(£ , 20 = 1) < ^dyn- Assuming that the particle es- 
cape time is comparable to the dynamical time, when the 
condition in equation [5] is met, it implies that the spec- 
trum of escaped particles is similar to that determined 
by the acceleration mechanism. We further discuss the 
implication of this condition in section IIIII below. Com- 
bining equations U and one obtains an upper limit on 
the value of the magnetic field at the acceleration site, 



B < AT 2 A 4 Z- 



-3 rpob 
^20 



V 1 G. 



(6) 



A lower limit on the value of the magnetic field can only 
be obtained once the acceleration radius r is specified. 
This will be discussed in section Hill below. 

In addition to synchrotron energy losses, energetic par- 
ticles can in principle lose their energy by interacting with 
the ambient photon field and with other nuclei. Inter- 
action with the photon field can result in energy losses 
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by Compton scattering (which results in negligible en- 
ergy losses for UHECR's, and therefore will not be con- 
sidered here), photopair production (Bethe-Heitler pro- 
cess), photo-production of mesons (mainly pions), and 
photodisintegration of the nuclei. Therefore, the condi- 
tion < acc < £ C ooi as given in equation[2]should be modified. 



Constraints from Disintegration of Energetic 
Nuclei 



The threshold energy for photopair production, e 



2m e (l + m e /m p ) ~ lMeV, is lower than the threshold 
energy for pion production, e£ h = m„(l + m„/2m p ) ~ 
145 MeV. The relative contributions of these two pro- 
cesses to the energy loss of a relativistic particle is deter- 
mined by the product of the cross section and inelasticity 
coefficient of the relevant process. For energetic protons, 
the product of cross section and inelasticity in pion pro- 
duction process is two orders of ma gnit ude larger than 
the product in pair creation process [66l . . Exact cal- 
culation must take into account the target photon spec- 
trum. It can be shown that for a photon number spec- 
tral index a — 2 photomeson production is somewhat 
more dominant, however both time scales are compara- 
ble (68|. For heavy nuclei this ratio drops by a factor 
Z 2 /A (see Refs. [HI, HH and references therein), which 
for iron nuclei still implies that photomeson production is 
more important than photopair production as an energy 
loss channel. 

However, the main energy loss channel of energetic par- 
ticles is photodisintegration [7(|. The threshold energy 
for this process is ef h ~ 10 MeV, larger than the threshold 
energy for photopair production. Here, we thus focus on 
photomeson production and photodisintegration as the 
main energy loss channels of energetic particles interact- 
ing with the photon field. 

An energetic nucleus having Lorentz factor 7^4 (in 
the plasma comoving frame) propagating through an 
isotropic photon background with differential number 
density n(£ 7 )efe 7 (at the energy range £ 7 ...£ 7 + <fe 7 ) loses 
energy by photodisintegration and photomeson produc- 
tion at a rate (see, e.g., Refs. [Til. Wi\) 



dis,7r 



n(e") 



de". 



'Ill A 



(7) 

Here, (Ja.tt{^') are the cross sections for photodisintegra- 
tion and photomeson production, and k(c') is the inelas- 
ticity coefficient in the photomeson production process. 
In calculating survival of one specie of heavy nuclei due 
to the photodisintegration process, equation [7] can still 
be used by inserting k(e') = 1. 

We estimate the differential photon number density 
at the acceleration radius r acc in the following way: as 
a rough approximation, the spectrum of RQ AGNs in- 
creases as vL v oc v 2 at low energies, v < v$ b = 10 13 Hz, 
and is flat (yL v oc u ) at higher frequencies, up to the 



X-rays « b ax ~ 10 18 Hz; See Refs. [H M, Iz3) & 
The differential photon number density can therefore be 
written as a broken power law, 



n(e) = n (e/e y 



(8) 



where a — for e < £0 an d a — 2 for e > Eq . Here, 
eo = hv^ b /r = 6.6 x 10 _14 r _1 erg. The normalization 
constant, no is found by normalizing to the total bolo- 
metric luminosity L£ b ol = 10 43 L^ 143 ergs -1 . Note that 
this bolometric luminosity does not correspond to the 
photon luminosity of the weak jets themselves, which 
are necessarily weaker. However, the survival of UHE 
nuclei depends on the entire photon field. Using the 
(comoving) energy density at the acceleration radius, 
u = L° b ol /(<inr 2 r 2 (3c) = J n{e)ede, one finds 



71-0 47regr 2 r 2 /3c[l/2+log(£ max /e )] 



6xl0 56 
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cm 3 erg 
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where we approximated log(e max /eo) ~ 10. 

For the photon spectrum given by equation^ the inner 
integral in equation [7] can be written as 



n(e") „ _ n /e ( 



(1 + a) V27A£o 



-(i+q) 



(10) 



(Note that the above result is accurate for e' j2 r y A > £q, 
while for e' /2ja < £0 there is a second term that can 
be neglected for large enough values of ja', see further 
discussion below.) 

In order to estimate the outer integral in equation [3 
we discriminate between photodisintegration and pho- 
tomeson production processes. In the photodisintegra- 
tion process, the main contribution to the outer integral 
in equation [7| is from photons in the energy bandwidth 
of the giant dipole resonance, whose energy is given by 
e GD R = 42.65 xi-°- 21 MeV for A > 4 [76]. Numerical fits 
to the experimental data gives the energy bandwidth of 
the resonance Aqdr = 8 MeV, and the maximum value 
of the cross section to be (Jq^a = 1A5A x 10 _27 cm 2 . Ap- 
proximating the outer integral by the contribution from 
the resonance, one finds 



cri <T ,AAGDR / 2-fASQ ' 
\ «GDR 



-(!+<*) 



7a(1+«) 



(11) 



Note that this result is in good agreement with the more 
detailed numerical calculation at ja < £gdr/2£q. (see Ref. 
[H| and figure [J) [97|. 

The photodisintegration time can be compared to the 
dynamical time, r/Tf3c [98], to obtain a lower limit on 
the acceleration radius which allows survival of the ac- 
celerated nuclei. The result depends on the target photon 
spectrum. For the assumed spectral index a ~ 2 above 
Vn b , the dissipation time obtains its minimum value at 
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J A ~ £gdr/2£o- For nuclei with this value of the Lorentz 
factor, one finds the survival condition, 



id 



t^i < 1 - r>2.1xl0 16 A 121 Lg b o , 43 < 13 



For 
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/3 cm. 
(12) 

nuclei, this gives r> 3 x 
~ ? r- 4 /3- 2 cm. Therefore, for 10 



parsec-scale weak jets (assuming here (3 ~ 0.3) of RQ 
AGNs, we can expect that heavy nuclei can not only be 
accelerated but also survive without significant photo- 
disintegration. We further note that for nuclei at higher 
energies, numerical calculations (e.g., Ref. [3|) indicate 
roughly comparable constraints on the acceleration ra- 
dius. As at these high energies (> £gdr/2£o = 10 19 eV) 
the simple analytical treatment may be insufficient given 
the contribution from non-GDR resonance, we carried a 
numerical calculation. We present in figure [1] a compari- 
son of the analytical approximation of the energy loss 
time to a more accurate numerical calculation, as well 
as the cooling time scales due to additional phenomena. 
While some deviation at the very high energies exist, 
clearly the analytical approximation used here is very 
accurate at lower energies, and does not deviate much 
from the exact numerical solution above 10 19 eV and 
below 10 20 eV, for the parameters values chosen. 

The second energy loss channel for energetic nucleons 
or nuclei is the photomeson production. The main con- 
tribution to the outer integral in equation [7] is from pho- 
tons at energies e pca k ~ 0.3 GeV, where the cross section 
peaks at the A-resonance, whose width is Ae w 0.2 GeV. 
For protons, the peak of the cross section at the A- 
resonance is <7 pea k ~ 5 x 10~ 28 cm 2 , and the inelasticity 
is tpcak ~ 0.2 [7l|. For heavier nuclei, the cross section 
is roughly proportional to A 17711 . while the inelasticity is 
proportional to A -1 (see Ref. (Zjl). Repeating the same 
calculation as for the photodisintegration process pre- 
sented above, one finds that the limitation on the accel- 
eration radius arises for particles at Lorentz factor 7^4 ~ 
e P cak/2£o- For these particles, we have t~ x ~ (2/(1 + 

a))n £ co'poakKpcak(Ae/epcak)(7A/0.5ep Ca keo 1 ) Q [II]. 
The requirement that td yn t^ 1 < 1 leads to a much looser 
constraint than that for photodisintegration. Note that 
this condition does not depend on A for a given ja ■ We 
thus conclude, that under conditions that enable survival 
of heavy nuclei, photomeson production is inefficient (see 
alsoRefs. 0,[6gj]). 

The number density of particles at this distance, n w 



L/4TTr 2 T 2 (3m p c 3 « 2x10" 



ob 

J bol,43 ' 18 



plies that spallation is not important as an energy loss 
channel of the energetic particles. The rate for spalla- 
tion process due to collision with other nuclei can be 
estimated as 



where a sp = 



t S p = cr sp r 

-26 



(13) 



5 x lO~'^ b A 2 / 3 cm 2 is the cross section for 
spallation of a nucleus (6f| ■ One can thus conclude that 
for acceleration at parsec scales, the spallation loss time 
is thus much longer than the dynamical time. 
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FIG. 1: The acceleration time scale and various cooling 
time scales of iron for RQ AGNs. Used parameters are 
Lti = 10 43 ergs" 1 , L B = e B L = 10 43 ergs' 1 , ug b = 10 13 Hz, 
photon spectral index a = for v < vo and a = 2 for 
v > vq, acceleration radius r = 10 19 cm, jet velocity (3 — 0.3, 
and r) = /3~ 2 = 10. Shown are the inverse of the relevant 
time scales. For the parameters shown here, the accelera- 
tion time (thin red line) is much shorter than the dynamical 
time (thick red line) and the synchrotron (dash-dotted red 
line) and Compton (dashed green line) cooling times of iron 
nuclei at energies below 10 20 eV. The thick, dash blue line 
shows the exact numerical results of the calculation of iron 
photodisintegration time scale. This calculation includes the 
nonresonant effect, and is similar to the ones presented in Ref. 
[l4| . The thin blue line, which by large overlaps the thick blue 
line, shows the analytical result discussed here, using equation 
1111 Note that the parameter set chosen allows for survival of 
ultrahigh-energy iron nuclei up to ~ 10 20 eV. 



III. ENERGETICS OF RQ AGNS AND 
EFFICIENCY OF UHECRS ACCELERATION 

In the previous section, we have concluded that the 
conditions inside the weak jets of RQ AGNs enable the 
acceleration and survival of UHE nuclei, if particle accel- 
eration occurs at parsec scale radii. Such radii may be 
characteristic of internal or standing shocks in a jet, or 
shocks produced by interaction of the jet and the ambient 
medium. In this section, we compare the energy budget 
in RQ AGNs with the energy requirement for UHECRs 
acceleration. 

The observed flux of UHECRs on earth is J w 5 x 
10~ 16 m~ 2 s _1 sr _1 (see Ref. |2J and references therein). 
The production rate of UHECRs within the GZK horizon 
of i? GZK » 100 Mpc [rmJ is therefore M = A-kJ/R gzk w 
1.7 x 10 36 Mpc -3 yr _1 . For energies above 10 19 5 eV, this 
implies energy production rate in energetic particles of 
E w 8.5 x 10 43 ergMpc~ 3 yr" 1 . This value is basically in 
agreement with more detailed estimates (see, e.g., Refs. 

PED- _ _ _ 

Following Refs. [43|, [57|, l78| (and references therein), 
we consider a typical RQ AGN jet luminosity L- Jct ~ 
10 42-43 erg s _1 . These values are based on measurements 
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in the optical band [78|; The typical radio emission is 
much weaker, ~ 10 38-39 ergs -1 , due to both synchrotron 
self absorption and free-free absorption. Although the 
value for Lj et is highly uncertain, we find this luminosity 
to be equal to or somewhat fainter than the magnetic 
luminosity required for acceleration of UHE nuclei (say, 
with Z > 10, see eq. The number density of RQ 
AGNs in the local universe is (^42 ~ 10 -3 Mpc -3 which 
is about 10 times larger than RL AGNs (see, e.g., Ref. 
[7{^]). Therefore, we expect that the energy input from 
RQ AGNs is ~ 3 x 10 46 ergMpc -3 yr" 1 . 

The number density of RQ AGNs drops with their lu- 
minosity [IS [M[. For example, the number density of 
AGNs with soft X-ray luminosity, Lx — 10 44 ergs _1 , is 
approximately ^44 ~ 10~ 6 Mpc~ 3 . While at present the 
jet luminosity is not easily deduced from the X-ray lu- 
minosity measurements (e.g., absorption has a different 
effect on measurements at the soft and hard X band, see 
Ref. 82]), we can still conclude that expected energy 
output from RQ AGNs jets significantly exceeds the en- 
ergy requirement for UHECRs sources. Therefore, RQ 
AGNs jets are energetic enough to supply acceleration of 
particles to ultrahigh energies. 

The estimated number density of RQ AGNs, $42 ~ 
10~ 3 Mpc~ 3 , implies that for isotropic distribution of 
sources, the rate of energy production in UHE particles 
in a single source should be roughly w 3 x 10 39 ergs _1 
at 10 19 ' 5 eV, in order for these objects to be the sources 
of UHECRs (and assuming that no additional sources of 
UHECRs exist). This value is less than a percent of the 
expected typical jet luminosity of these objects. 

This obtained value of ~ 10 ergs -1 , while well 
within the energy limit of a typical RQ AGN, should 
be considered as a lower limit. First, the estimated num- 
ber density of RQ AGNs does not take into account the 
possibility that in some of these objects the jet luminosi- 
ties may be insufficient (see further discussion below). 
Second, the uncertainty that exists in the value of the 
efficiency parameter r), as well as in the shock velocity (3 
implies, via equation [2] that the required luminosity for 
acceleration of iron nuclei to 10 20 eV may be somewhat 
higher than the value 10 42 ergs -1 considered here, in 
which case the number density of available sources may 
be lower. Therefore, even though our scenario remains 
essentially viable, one should keep in mind the high un- 
certainty that exists in the details of the physical pro- 
cesses. Nonetheless, the values obtained here can give a 
first approximation of the physical requirement from the 
source of UHECRs. The fact that they are well within 
the known energy budget of these objects leaves an am- 
ple margin to incorporate many of the uncertainties in 
the values of the physical parameters. 

The above estimate does not take into account accel- 
eration of particles to lower energies, which accompa- 
nies the production of UHECRs. Therefore, the required 
cosmic-ray input per source depends on the spectral in- 
dex of the accelerated particles in the source. Although 
very steep source spectra can lead to an 'energy crisis' 



(see, e.g., Refs. (HHHI), RQ AGNs can be viable sources 
if sufficiently flat spectra are achieved. As both the jet 
luminosity and number density are uncertain at present, 
the results obtained here can perhaps serve as a guide- 
line. Additional source of uncertainty lies in the energy 
losses caused by shock waves formed if a large number of 
charged particles escape from the source. While the en- 
ergy of the escaping particles is reduced, still the highest 
energy particles, for which the escape time is the short- 
est, will be less affected. The exact details of this effect 
depend on the exact spectrum of the escaping particles, 
and are thus beyond the scope of this manuscript. How- 
ever, we point here that if only the highest energy par- 
ticles escape, this effect is not expected to play a signifi- 
cant role. Then, the overall cosmic-ray energy spectrum 
may be achieved as a superposition of contributions from 
many RQ AGNs with different cosmic-ray luminosities 
and maximum energies. Finally, our model requires a 
significant fraction of heavy element composition in the 
jet. As the jet material is ejected from the inner part of 
the disk, existence of metals in the disk inevitable leads 
to their existence in the jet, unless nuclei are disrupted 
in the jet base and the accretion flow in the disk. 

The constraints found in EjTT] above imply that UHECR 
production is possible if the acceleration radius is of the 
order of parsec, in order to avoid energy loss by photo- 
disintegration. We note that the calculations carried in 
[jlT] considered the bolometric luminosity, which is higher 
than the jet luminosity. As a result, the constraints de- 
rived above are more restrictive than the constraints that 
would have obtained by considering the (uncertain) lu- 
minosity in the jet only. 

Using the results derived in $H] we can obtain a 
lower limit on the value of the magnetic field at the 
acceleration site if we assume the acceleration radius 
r, and characteristic velocity, /3. If UHECRs are 
composed of iron nuclei, assuming (3 = 0.3, equa- 
tion [TJ] requires acceleration radius of r ~ 10 pc. 
At this radius, the Hillas condition (equation [1]) im- 

plies B>10- 4 77 /3 3 1/2 ^2 oiboi,43<i3~ 2 G. For Car- 
bon nuclei, the lightest nuclei that are consistent with 
the restrictions on the luminosity obtained in equa- 
tion^ the minimum acceleration radius consistent with 
equation [T3] should be at 5 x 10 18 cm. The mag- 
netic field at this radius should thus exceed B > 3 x 
10~ 3 rj (3 3 _ 1/2 E°l L° b olA3 vg 13 ~ 2 G. We note that for both 
scenarios, sub-Gauss values of the magnetic field is con- 
sistent with the requirement that the cooling time is 
shorter than the dynamical time (see equation [5]), in 
which case the spectrum of the escaping particles is sim- 
ilar to that determined by the acceleration mechanism. 
We can thus constrain the required value of the magnetic 
field at the acceleration site to be in the range 10 -3 — 1 G. 

This ~ 1 parsec acceleration radius derived here is con- 
sistent with the observed scale of weak jets observed in 
nearby Seyfert galaxies HI, US US USUI]- At this radius, 
internal interactions with the jets, or interaction of the 
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jet with the interstellar medium material are expected to 
produce shock waves, that dissipate the jet kinetic en- 
ergy. These shock waves are plausible acceleration sites 
for energetic particles. The lower limits on the values of 
the magnetic field derived above can be translated into 
lower limits on the equipartition value of the magnetic 
field at the acceleration site: using B 2 /8ir = esu (see 
discussion above equation [2~]), one finds that for iron nu- 
clei, eg > 4 x 10 -3 , while for carbon nuclei, €b > 10 -1 . 
These values are both below equipartition, and are con- 
sistent with values of the magnetic field expected to be 
produced in shock waves. 



IV. OBSERVATIONAL IMPLICATIONS 

In the previous sections, it has been shown that RQ 
AGNs can be viable sources of UHECRs, provided that 
the UHECRs composition is dominated by heavier nu- 
clei (e.g., iron) rather than protons. This claim is based 
on two general considerations: (1) For a heavy nuclei 
composition, the required magnetic luminosity is weak 
enough to enable RQ AGNs to accelerate UHECRs; and 
(2) RQ AGNs are about ten times more numerous than 
RL AGNs, so that they can significantly contribute to 
the observed UHECR flux. 

If the assumption presented here is correct, in addi- 
tion to the contribution of RQ AGNs to the observed 
UHECR flux, we expect a contribution from RL AGNs. 
This is because the same acceleration mechanism can be 
expected to work in the latter. Although it is difficult to 
estimate the relative contribution of these two classes of 
objects, we expect that heavy nuclei (presumably, iron) 
from RQ AGNs can mask the correlation that may be 
found between the arrival direction of UHECRs and RL 
AGNs. In this respect, our scenario is consistent with the 
recent "disappointing model" [12] , which does not predict 
a significant correlation of the most energetic CR's with 
AGNs. 

Nonetheless, the conditions within the jets of RL AGNs 
may enable proton acceleration to ultra-high energies 
within these objects as well. Thus, if our scenario is 
correct, even if the correlation that can be expected be- 
tween the arrival direction of UHECRs and positions of 
RL AGNs is reduced by the additional contribution from 
RQ AGNs, some correlation between the arrival direc- 
tions of UHE protons and the positions of RL AGNs is 
still expected. The strength of this correlation depends 
on the efficiency of particle acceleration in RQ and RL 
AGNs, and on the composition of particle acceleration 
in RL AGNs: as our model allows, in principle, proton 
acceleration to UHE in RL AGNs, the observed UHE- 
CRs flux may be composed of two distinctive popula- 
tions: heavy nuclei accelerated in RQ AGNs, and pro- 
tons accelerated in RL AGNs. A testable consequence of 
this idea is that UHECRs whose arrival directions cor- 
relate with the position of RL AGNs will show, on the 
average, lighter composition than UHECRs whose arrival 



direction is not correlated with RL AGNs. Nonetheless, 
we stress that this is highly uncertain, due to two main 
reasons: first, RL AGNs may accelerate heavy nuclei to 
high energies as well; and second, the efficiencies of parti- 
cle acceleration and escape from these sources are highly 
uncertain. 

We can therefore conclude that our model does not 
rule out some correlation between the arrival direction of 
UHECRs and positions of RL AGNs. We note though, 
that the current observational status is highly uncertain. 
While some correlation was reported [3l|, |37( , its strength 
seems not to be fully determined yet (the recent results 
reported by Ref. [37 1 in dicate much weaker correlation 
than earlier reports 31(). Therefore, future anisotropy 
search focusing on UHE proton events should be very 
important in confirming this idea. 

In our scenario, the detection of secondary gamma rays 
and neutrinos from individual RQ AGNs whose luminos- 
ity is thought to be L ~ 10 42- ergs -1 , seems to be 
difficult. This is due to the fact that, as we saw in W A[ 
in the weak jets in RQ AGNs photodisintegration is the 
dominant energy loss channel, hence copious production 
of energetic 7r's is not expected. We leave the details of 
this calculation for future work. 

On the other hand, the typical luminosity of jets of RL 
AGNs can be a few orders of magnitude brighter than 
RQ AGNs. Therefore, RL AGNs are more favorable tar- 
gets for the purpose of detecting energetic secondaries. 
In particular, secondaries originating from very powerful 
AGNs like FR II galaxies (see, e.g., Ref. .84J and refer- 
ences therein) or very nearby RL AGNs such as Cen A 
(e.g., Ref. [23| and references therein) might be detected 
by Fermi or IACTs. Clusters of galaxies hosting power- 
ful AGNs might also be viable candidate for secondary 
detection surveys (see, e.g., Refs. [iHHH and references 
therein). 



V. SUMMARY AND DISCUSSIONS 

In this paper we have considered radio quiet AGN's 
as possible sources of UHECRs. So far, jets in these ob- 
jects were not considered as sources of the highest energy 
cosmic rays around 10 20 eV, since they are not luminous 
enough to support acceleration of protons to these ener- 
gies. However, the recent results of the AUGER collabo- 
ration shows indications that at the highest energies the 
composition of cosmic rays may be dominated by heavy 
nuclei [H, . The assumption that UHECRs are heavy 
nuclei eases the constraint on the source luminosity (see 
eq. [2]). Thus it allows, in principle, radio-quiet AGNs to 
be the sources of UHECRs. 

We have calculated in JjTT] the constraints on the ac- 
celeration site which enable the acceleration of heavy 
nuclei to high energies. The most restrictive constraint 
arises from the photo-disintegration process. We showed 
in ^11 At that for typical nearby AGN's with bolometric 
luminosity L^ ol « 10 43 ergs -1 , energetic nucleus can sur- 
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vive photo-disintegration if the acceleration takes place 
on a parsec scale (see eqs. [T^]). Interestingly, this is the 
same scale at which weak jets were seen in these objects 
[541 [55l . HH, [58|. This fact further supports our idea, 
since interaction of jets with the surrounding material 
inevitably leads to creation of shock waves, which are 
the most plausible acceleration site of particles to high 
energies. 

The question of the acceleration of particles to ultra- 
high energies may depend on the shock velocity. In the 
Bohm approximation r\ oc /3 -2 , some assumptions about 
the shock velocities in the jets are needed. For example, 
a shock velocity of (3 ~ 0.1 requires magnetic luminos- 
ity of €bL> 10 43 ' 5 ergs -1 in order to enable acceleration 
of iron nuclei to the maximum observed energy (see CHI 
equation^, which may be too low. However, somewhat 
higher values, (3 ~ 0.3 may be sufficient. Since currently 
there is a high uncertainty in the determination of the 
shock velocities in RQ AGNs, it is not possible to put 
further constraints. We have also showed in ^Illl that the 
minimum inferred values of the magnetic field required 
to confine the accelerated particles at parsec scale jets is 
sub-Gauss. This value was shown to be smaller than the 
equipartition value, which is also consistent with mod- 
els of particle acceleration in shock waves that may be 
produced on this scale. 

We further showed in i jllll that the required luminosity 
in UHECRs could be less than a percent of the total bolo- 
metric luminosity of more abundant nearby AGNs. If the 
spectrum of produced cosmic rays is a power law with in- 
dex close to d log N/d log E sa —2, as suggested by mod- 
els of particle acceleration in non-relativistic shock waves 
(63| . then the total luminosity in energetic particles (at 
all energies) could still be much smaller than the total 
bolomctric luminosity. Even for higher power law index, 
p = 2.3, we find that the total energy requirement for ac- 
celeration of cosmic rays above GeV is « (10 10 ) 3 = 10 3 
times higher than the energetic requirement from UHE- 
CRs alone, which is still (marginally) consistent with the 
total energy budget in RQ AGNs (see discussion in EjHII) . 
Assuming that the main radiative source in these objects 
is the accretion disk, this implies a high efficiency, of the 
order of tens of percents, in the conversion of accretion 
energy to acceleration of particles. Since large uncertain- 
ties exist in the efficiencies of both the energy conversion 
in the jet and the acceleration of particles, we can only 
conclude that our model is consistent with the energetic 
requirements, provided a high efficiency is achieved in 
both these processes. 

Spectral synthesis as well as chemical enrichment mod- 
els predict that AGN's are metal rich. The metallicitics 
in the broad emission line reg ion is typically ~ 1 to > 10 
times the solar metallicity [8(| [87} , and grows with the lu- 
minosity (88| . Thus, abundant existence of heavy nuclei 
in AGN's disks is expected. As the plasma jet is com- 
posed of material from the disk, abundant population of 
heavy nuclei in the jet is expected. Thus, AGN's may be 
one of the natural sources of high energy heavy nuclei. 



If indeed the acceleration takes place on a parsec scale, 
then the main energy loss channel is photodisintegration, 
rather than photomeson production. As discussed in mVl 
this fact implies that we do not expect abundant pro- 
duction of energetic neutrinos, that may result from the 
decay of energetic pions. On the contrary: the results 
in equations [T^] indicate that the relative contribution of 
photomeson production may be no more than few per- 
cents of photodisintegration as an energy loss channel of 
UHECRs. Thus, we do not expect a copious production 
of neutrinos under the conditions which enable the ac- 
celeration and survival of heavy nuclei, as considered in 
this paper. 

After they escape from their sources, UHECR nuclei 
can be subject to photodisintegration in the intergalactic 
space. UHE nuclei with energy > 10 19 eV mainly inter- 
act with the cosmic infrared background (CIB) photons 
[ESIHlIsi- UHE iron nuclei with energy < 10 20 eV have 
a mean free path of > 500 Mpc, while the mean free path 
of oxygen nuclei at a similar energy is only ~ 30 Mpc 
|9Cj |. Both increase rapidly as the energy of the nuclei 
decreases. Thus, the GZK horizon of UHE nuclei is com- 
parable to, or even somewhat larger than the GZK hori- 
zon of energetic protons (see also Ref. [l|). 

To summarize, we have pointed out here that the as- 
sumption that UHECRs are composed of heavy nuclei, as 
is suggested by recent AUGER results [H, [33] , enables 
radio quiet AGN's to be their sources. This picture is 
consistent with what is currently known about the exis- 
tence of parsec scale jets seen in some of these objects, 
as well as the energy constraints. Moreover, this pic- 
ture is supported by the recent analysis carried by Ref. 
[oil ] , in which a strong correlation between the arrival di- 
rections of UHECRs and the sky coordinates of AGN's 
detected by the Swift-BAT within the GZK cutoff (which 
are largely radio-quiet) was found. As more data are col- 
lected by the AUGER collaboration, a clearer picture of 
the composition of UHECR's will become available, al- 
lowing further constraints on possible sources. 
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